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The TeV blazar Markarian 421 underwent multi-TeV flaring during April 2004 and simultaneously 
observed in the X-ray and TeV energies. It was observed that the TeV outbursts had no counterparts 
in the lower energy range, which implies that this might be an orphan flare. We show that Fermi- 
accelerated protons of energy < 168 TeV will interact with the low energy tail of the background 
synchrotron self-Compton photons in the inner region of the blazar to produce the multi-TeV flare 
and our results fit very well the observed spectrum. Based on our study, we predict that the blazars 
with a deep valley in between the end of the synchrotron spectrum and the beginning of the SSC 
spectrum are possible candidates for orphan flaring. Future possible candidates for this scanario are 
the HBLs Mrk 501 and PG 1553+113 objects. 


I. INTRODUCTION 

Active galactic nuclei (AGN) emit electromagnetic ra¬ 
diation from radio to gamma-rays and exhibit large lu¬ 
minosity variations on time scales ranging from less than 
an hour up to several years. A super massive black hole 
is believed to sit at the center of the AGN surrounded 
by an accretion disk in the inner region and a torus of 
gas cloud in the outer region. Oppositely directed rel¬ 
ativistic jets are ejected from the AGN which are per¬ 
pendicular to the accretion disk and the torus. In the 
framework of the unification scheme of AGN, blazars and 
radio galaxies are intrinsically the same objects, viewed 
at different angles with respect to the jet axis. When 
the angle between the jet and the line of sight is small 
it is called blazar and in contrast, for radio galaxies, the 
angle between the jet and the light of sight is large. Al¬ 
most all AGN detected at very high energy (VHE) (> 100 
GeV) are blazars with the exception of the three objects, 
Centaurus A (Cen A) [ill, 0, M87 and NGC 1275 which 
are radio galaxies [3, b|. The spectral energy distribution 
(SED) of these AGN have a double peak structure in 
the v — vF v plane. The low energy peak corresponds to 
the synchrotron radiation from a population of relativis¬ 
tic electrons in the jet. Although the general consensus is 
that the high energy peak corresponds to the synchrotron 
self Compton (SSC) scattering of the high energy elec¬ 
trons with their self-produced synchrotron photons, this 
result remains inconclusive for various reasons. However, 
the so called leptonic model is very successful in explain¬ 
ing the multi wavelength emission from blazars and FR 

I galaxies [l 45}. 

Although the SSC scenario seems to work very well to 
explain the SED of AGN up to the second peak , 
difficulties arise in explaining the multi-TeV emission de¬ 
tected in Cen A0, flares from the radio galaxy M87 
the flares from blazars 1ES 1959+650 0, [7o| and Markar¬ 
ian 421 (Mrk 421) [ll| . Also the inevitable outcome of the 
leptonic scenario is that, emission in multi-TeV energy 
has to be accompanied by a simultaneously enhanced 


emission in the synchrotron peak. Unfortunately the en¬ 
hanced synchrotron emission was not observed in the flar¬ 
ing of 1ES 1959+650 in June 2002@ and also probably 
in the flaring of Mrk 421 in April 2004[ll|, which im¬ 
plies that the SSC model may not be efficient enough to 
contribute in the multi-TeV regime. 

To explain the orphan flaring of 1ES1959+650 a 
hadronic synchrotron mirror model was proposed by 
Bottcher [12| • In this model, the high energy protons from 
the jet interact with the primary synchrotron photons 
that have been reflected off clouds located at a few pc 
above the accretion disk. These photons are blue shifted 
in the jet frame so that there will be a substantial de¬ 
crease in the high energy protons to overcome the thresh¬ 
old for A-resonance. Similarly a structured leptonicjet 
model is proposed to explain the orphan TeV flare fl3j. 

In the hadronic models, the second bump is believed 
to be formed due to the synchrotron emission from the 
ultra-high energy protons, or emission from the sec¬ 
ondary electron-positron pairs or muons produced from 
the charged pion decay as a reasult of the interaction 
of the high energy protons with the background low en¬ 
ergy photons. In this process the decay of neutral pions 
to 7-ray pairs can explain the high energy peak in the 
multi-TeV range. 

Previously it was shown that the above hadronic pro¬ 
cesses are inefficient and to explain the high energy peaks 
efficient accleration of relativistic protons to ultra-high 
energies within the inner part of the jet outflow is re¬ 
quired. At the same time the jet kinetic power has to 
exceed the Eddington luminosity by orders of magnitude. 
The hadronic model is employed by Zdziarski et al.flfj 
to explain the broad band spectra of radio-loud active 
galactic nuclei and by Cao et al.[l5j to explain the TeV 
spectrum of the blazar 1ES 1101-232. In both these cases 
it is shown that super Eddington luminosity in proton is 
required to explain the high enery peaks. This situation 
arises because the photon density in the jet is low. How¬ 
ever, this problem can be circumvented if one assumes 
that within the jet there is a dense inner jet region [l6|. 
The multi-TeV peak can also be produced by pp interac- 
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tions. It was shown earlier that normally the pp process 
is inefficient to produce 7 -rays in the blazar environment 
unless the jet-cloud interaction is taken into account [l7j- 
lisij l . The emisson of gamma-rays by the interaction of a 
population of massive stars surrounding the AGN jets is 
studied in[20L f2li|. Also emission of high energy gamma- 
rays by the interaction of the AGN jets with the tidally 
disrupted atmosphere of Red Giant is studied i 11 fill. [231 . 


II. PHOTOHADRONIC MODEL 


To address the orphan TeV flaring of 1ES 1959+650 
and the multi-TeV emission from Cen A and M87, Sahu 
et al. [HIHIIHH have used the hadronic model. Recently 
this model is also used to estimate the neutrinos flux 
from the TeV-blazars and their spatial correlation with 
the IceCube eventsf26|. Here we use the same hadronic 
model to explain the multi-TeV orphan flaring of Mrk 
421. In this model the Fermi-accelerated high energy 
protons interact with the SSC photons in the core region 
of the jet (few times the Schwarzscliild radius Rs) to 
produce the A-resonance. Subsequently the A-resonance 
decays to charged and neutral pions as follows, 


p + 7 —>- A + —>■ 


p 7T W 


(1) 


The decay of neutral pions to TeV photons gives the 
multi-TeV SED. Throughout our work we use natural 
units c = h = 1. The relationship between the 7 r°-decay 
TeV photon energy E 1 and the target SSC photon energy 
e 7 in the observer frame is given by 

V 2 

£Le 7 ~ 0.032 --GeV 2 , (2) 

77 (1 + z) 2 

where V is the Doppler factor and 2 is the redshift. The 
observed TeV 7 -ray energy and the proton energy E p are 
related through 


Ep — p Ery , (3) 

where T is the bulk Lorentz factor of the relativistic jet. 

Multi-TeV emission is observed from many blazars and 
FRI galaxies by Cherenkov telescope arrays. These emis¬ 
sion are classified into non-flaring and flaring events. For 
the non-flaring events (e.g. emission from Cen A0), 
the injected proton spectrum is a power-law given by 
dN p /dE p oc E~ a with a > 2. These protons will inter¬ 
act with the background SSC photons having comoving 
number density ra 7 (henceforth ' implies the jet comov¬ 
ing frame) and satisfying the kinematical conditions in 
Eqs.© and ([3]). As discussed in Ref, fit} . the flaring oc¬ 
curs within a compact and confined volume of a smaller 
cone which is enclosed in a bigger cone as shown in Fig. 
H] The internal compact jet models are also proposed to 
explain the fast variability of the blazars Mrk 501 and 



FIG. 1. Geometry of the orphan flaring of blazar Mrk 421: 
the hidden interior compact cone (jet) is responsible for the 
multi-TeV orphan flaring and the exterior cone corresponds 
to the normal jet. 


PKS 2155-304(273, [H|. In our case, the injected proton 
spectrum is a power-law with an exponential decay. In a 
unified manner we can express the injected proto n sp ec- 
trum for both non-flaring and flaring processes as|16l 1291] 


dNp f 1 , non-flaring 

dE p ^ p \ e~ Ep / Ep - c , flaring 


(4) 


where E PtC is the break energy for the high energy pro¬ 
tons and above this break energy the proton spectrum 
falls very fast. The high energy protons will interact in 
the flaring region having the comoving photon number 
density n' ^ to produce the A-resonance. The photon 
density in the flaring region is much higher than the rest 
of the blob (non-flaring region with the photon density 
n 7 ) probably due to the copious annihilation of electron 
positron pairs, splitting of photons in the magnetic field, 
enhanced IC photons in this region and Poynting flux 
dominated flow which can form from the magnetic re¬ 
connection in the strongly magnetized plasma around the 
base of the jet [27], 130, Dll ■ Here we assume that the spec¬ 
trum has a exponential decay above the E PtC only for the 
flaring case. In this model the inner compact region has 
a size R’j smaller than the size of the outer region R' b i.e. 
R'f < R' b . Unfortunately, we do not know exactly the 
properties of the inner jet except for the number den¬ 
sity, the energy density of photons and and the magnetic 
field are higher than the corresponding values in the outer 
jet. Due to higher photon density in the inner region, the 
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multi-TeV 7 -rays and neutrinos can be produced through 
intermediate A-resonance as shown in Eq. (fljl. The op¬ 
tical depth of the A-resonance process in the inner jet 
region is given by 


T vl = v! lf (j\B! s . (5) 

The efficiency of the py process depends on the physical 
conditions of the interaction region, such as the size, dis¬ 
tance from the base of the jet, photon density and their 
distribution in the region. 

In the inner region we compare the dynamical time 
scale t' d = Rj with the py interaction time scale t' pl = 
(n(. to constraint the seed photon density so 

that multi-TeV photons can be produced. For a mod¬ 
erate efficiency of this process, we can assume t' pl > t' d 
and this gives r P7 < 2 , where we take the inelasticity 
parameter K pl = 0.5. Also by assuming the Eddington 
luminosity is equally shared by the jet and the counter 
jet, the luminosity within the inner region for a seed pho¬ 
ton energy e 7 will satisfy (4™.^ jR'j e 7 ) -C Lsdd/ 2. This 
puts an upper limit on the seed photon density as 


n 


7> / 




F Edd 

8ttA/ e 7 


( 6 ) 


photon energies E 1 1 and E 1 2 can be expressed as 

= n 7( £ 7i) (^n\ + e -( E ^-E y2 )/E c /g\ 

Fj(E 12 ) n 7 (e 72 ) \E J2 J 

where Ey 1 correspond to the proton energy E P12 . In 
this derivation we have used the relations E P 1 /E P2 = 
Ey 1 /E l2 , and E pl 2 /E PtC = E 112 /E c . By taking the ra¬ 
tio we get rid of the proportionality constant. By using 
the known flux at a particular energy in the flaring/non¬ 
flaring state, we can calculate the flux at other energies 
using Eq. © and also the normalization constant can be 
calculated. 

In terms of SSC photon energy and its luminosity, the 
photon number density n' is expressed as 


n 7 (e 7 ) = 77 


L~f,ssc{ 1 + z) 
■D 2 +k 4t jR'lt-i ’ 


( 10 ) 


where p ~ 1 is the efficiency of SSC process and k ~ 
(0 — 1 ) depending on whether the jet is continuous or 
discrete. In this work we consider n = 0. The SSC 
photon luminosity is expressed in terms of observed flux 
(^ssc^) = e 2 dV 7 /de 7 ) which is given by 


From Eq.© we can estimate the photon density in this 
region and this is discussed in Sec. 3. On the other hand, 
the photon density in the outer region can be calculated 
from the observed flux. For this we assume a scaling 
behavior of the photon densities in the flaring and the 
non-flaring regions which we take to be the following,, 

</( £ 7i) = <(4) m 

n 7,/( £ 72) ^7 ( £ 72 ) 7 

this implies that, the ratio of photon densities at two 
different background energies e' and e' in flaring and 
non-flaring states remains the same. For a self consistent 
treatment, in principle we should use the photon den¬ 
sity n'^ j in the hidden internal jet and solve the coupled 
transport equations for leptons and photons along the 
jet axis by taking into account their respective cooling 
mechanisms as well as the injection spectrum of the pri¬ 
mary particles [32]. To avoid this complication we assume 
the scaling behavior of the photon densities in different 
background energies as shown in Eq. 0 . In the pho- 
tohadronic scenario, the number of 7r°-decay photons at 
a given energy are proportional to both the number of 
high energy protons and the density of the SSC back¬ 
ground photons in the jet, i.e. N(Ey) oc N(E p )n:' . For 
the flaring case n'y is replaced by the photon density in 
the flaring region given by n! 1 The 7 -ray flux from the 
7T° decay is then given by 


Ly,ssc 


4t rd|$ssc(e 7 ) 

( 1+*) 2 


( 11 ) 


Furthermore, by using Eq.©, we can simplify the ratio 
of photon densities given in Eq. 0 to 


_ 4 > gsc(£7i) E 1 1 
n '-f( e 72 ) 4 >ssc(e 72 ) E i2 


Now we can express Eq. © in terms of observed SSC flux 
and E 7 as, 


F i( E n) = ^ssc(e 7 i) /EyA 0+3 e -( Eyi -Ey 2 )/E c 
Fy(E 12 ) 4>ssc , (e 72 ) \Ey 2 J 

(13) 

Although Eqs.© and (fl3l) give the same result, the latter 
form is much more simpler than the former one. The one 
in Eq. (fl3l) uses the SED of SSC photon calculated using 
the leptonic model. Here the multi-TeV flux is propor¬ 
tional to E~ a+3 and 4 ) 33 ( 7 (e 7 ) while in the former case 
it is proportional to E~ a+ 2 and to the photon number 
density n 7 (e 7 ). Finally we would like to point out that 
in the photohadronic process (py), the multi-TeV photon 
flux can be expressed as 


F(Ey) = A 7 4 >ssc(e 7 ) 


( E i 

V TeV 


—a+3 


(14) 


Fy(Ey) = E, 


_ „ 2 dN{Ey) 


dE~ 


oc E, 


2 dN(E p ) ^ 


dE n 


7,/' 


( 8 ) 


Using the scaling behavior of Eq. (0, the observed multi- 
TeV photon flux from 7 r°-decay at two different observed 


where A 7 is a dimensionless constant and both e 7 and 
Ey satisfy the condition given in Eq.©. This formula 
will be used to calculate the multi-TeV flux from both 
non-flaring (without exponential decay term) and flaring 
events from AGN and their subclass if the emission is 
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FIG. 2. The SED of Mrk 421 is shown in all the energy bands 
which are taken from Ref. [Il] . The flare of April 2004 in 
multi-TeV energy is shown here. The hadronic model fit to 
the April 2004 data is shown as continuous line to the extreme 
right. The shaded region is the energy range of SSC photons 
where the Fermi-accelerated protons are collided to produce 
the A-resonance. 


due to photohadronic process from the core region. For 
different blazars/AGN, the value of A 7 will be different 
which we shall discuss in the next section. We can cal¬ 
culate the Fermi accelerated high energy proton flux F p 
from the TeV 7 -ray flux through the relation 

F P (E P ) = 7.5 x (15) 

Tp-y[F p ) 

The optical depth t p1 is given in Eq. (USD- For the ob¬ 
served highest energy proton energy E p , F p (E p ) will be 
smaller than the Eddington flux Fem ■ This condition 
puts a lower limit on the optical depth of the process 
and is given by 


Tp 7 (Ep) > 7.5 x 


F 1 {E 1 ) 

FEdd 


(16) 


From the comparison of different times scales and from 
Ea. (fTT)l) we will be able to constraint the seed photon 
density in the inner jet region. 


III. MRK 421 


Mrk 421 is a high synchrotron peaked BL Lac object 
(HBL) and is the first extragalactic source with a red- 
shift of z=0.031 to be established as a TeV emitter[33j]. 



Log E y [TeV] 

FIG. 3. The continuous curve is the hadronic model fit to the 
multi-TeV flaring data of Mrk 421. 


It has a luminosity distance d,L of about 129.8 Mpc. 
Its central supermassive black hole is asumed to have 
a mass Mbh — (2 — 9) x 10 8 M 0 corresponding to a 
Schwarzschild radius of (0.6 — 2.7) x 10 14 cm and the 
Eddington luminosity Lsdd = (2.5 — 11.3) x 10 46 erg s _1 . 
The synchrotron peak of its SED is in the soft to medium 
X-ray range and the SSC peak is in the GeV range. It 
is one of the fastest varying 7 -ray sources. In the past 
through dedicated multi wavelength observations, the 
source has been studied intensively. These studies show a 
correlation between X-rays and very high energy (VHE) 
7 -rays. A one-zone SSC model explains the observed 
SED reasonably welljUj. Several major flares were ob¬ 
served in the 2003/2004 season. During April 2004, a 
large flare took place both in the X-rays and the TeV 
energy band. The flare lasted for more than two weeks 
(from MJD 53,104 to roughly MJD 53,120). But due to 
a large data gap between MJD 53,093 and 53,104, it is 
difficult to exactly quantify the duration. The source was 
observed simultaneously at TeV energies with the Whip¬ 
ple 10 m telescope and at X-ray energies with the Rossi 
X-ray Timing Explorer (RXTE). It was also observed si¬ 
multaneously in radio and optical wavelengths. During 
the flaring it was observed that, the TeV flares had no co¬ 
incident counterparts at longer wavelengths. Also it was 
observed that the X-ray flux reached its peak 1.5 days 
before the TeV flux did during this outburst. So it is 
believed that, the TeV flare might not be a true orphan 
flare like the one observed in 1ES 1959+650. On the 
other hand remarkable similarities between the orphan 
TeV flare in 1ES 1959+650 and Mrk 421 were observed, 
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TABLE I. These parameters (up to B') are taken from the 
one-zone synchrotron model of ref. El which are used to fit 
the SED of Mrk 421. The last three parameters are obtained 
from the best fit to the observed flare data in our model. 


Parameter 

Description 

Value 

Mbh 

Black hole mass 

(2 - 9) x 10 s M q 

z 

Redshift 

0.031 

r 

Bulk Lorentz Factor 

14 

V 

Doppler Factor 

14 

K 

Blob Radius 

0.7 x 10 16 cm 

B' 

Magnetic Field 

0.26 G 

R'f 

Inner blob Radius 

3 x 10 il:> cm 

a 

Spectral index 

2.7 

E c 

7 -ray Cut-off Energy 

6.2 TeV 


including similar variation patterns in X-rays. 

By using the one-zone SSC model the average SED of 
Mrk 421 is fitted in Fig. 11 of ref. pjl- I n Fig. [2]we have 
shown the SED of Mrk 421. In this figure, the red crosses 
to the extreme left are the measurements in the radio 
frequency band measured by 26 m telescope at the Uni¬ 
versity of Michigan Radio Astronomy Observatory (UM- 
RAO) and by the 13.7 nr Metsahovi radio telescope at 
Helsinki. The single cross (red) in the optical range is 
measured by Fred Lawrence Observatory (FLWO) 1.2 nr 
telescope. The points of the first peak in the X-ray range 
are from Rossi X-ray Explorer (RXTE) El ■ 

In the one-zone leptonic model, the blob of size R' b ~ 
0.7 x 10 16 cm moves down the conical jet with a Lorentz 
factor T ~ 14 and a Doppler factor of T> = 14. The 
emitting region is filled with an isotropic electron popu¬ 
lation and a randomly oriented magnetic field B' = 0.26 
G. In the present work we study the flaring of Mrk 421 
during April 2004. We use the parameters of the one- 
zone leptonic model of ref. 11]. The parameters of the 
one-zone synchrotron model are summarized in Table HI 
In principle the Lorentz factor in the inner jet should 
be larger than the outer jet. But here we assume that 
Tout — Tin ~ T. 


IV. RESULTS 


The flaring of Mrk 421 in April 2004 was observed 
in the energy range 0.25TeU(6.0 x 10 25 Hz) < E 1 < 
16.85 TeU(4.1 x 10 27 Hz) by the Whipple telescope. 
In the hadronic model discussed above, this corre¬ 
sponds to the Fermi accelerated proton energy in the 
range 2.5 TeV < E p < 168 TeV and the correspond¬ 
ing background photon energy will lie in the range 
23.6MeV(5.7 x 10 21 Hz) > e 7 > 0.35MeU(8.4 x 
10 19 Hz). This range of e 7 is the shaded region shown in 
Fig. [5] which is in the low energy tail of the SSC photons. 
For the calculation of normalized multi-TeV flux we take 
into account one of the observed TeV fluxes from the 
flare with its corresponding energy and with the use of 
Ea. (fT51) calculate other TeV fluxes. In this model the 


free parameters are the spectral index a and the TeV 7 - 
ray cut-off energy E c which are adjusted to obtain the 
best fit and the values are a = 2.7 and E c = 6.2 TeV. 
This value of E c corresponds to the proton cut-off en¬ 
ergy E PjC = 62 TeV and the background SSC photon 
enery £y : ssc = 0.96 MeV (2.3 x 10 20 Hz) which is very 
close to the beginning of the SSC energy as shown in Fig. 
[2] This shows that for orphan flaring, the cut-off energy 
E c is due to the change from the synchrotron band to 
the SSC band and can be calculated from their crossover 
energy. 

In Fig. [3] we show both the observed multi-TeV SED 
and our model fit to it. In our results, the presence of 
^ssc(^) in Eq. (ThU) modifies the power-law with the ex¬ 
ponential fall scenario. From the best fit parameters we 
obtain the value of the dimensionless constant A 1 ~ 20 
in Eq. CE3D- By using the best fit parameters, we have 
also calculated the value of A 1 from the multi-TeV flare of 
1ES 1959+650 and M87 and the multi-TeV emission from 
Centaurus A, which are given as 86 , 1.86 and 8.6 x 10 -4 
respectively. We observe that for orphan flaring the con¬ 
dition A 7 1 is satisfied as is seen from 1ES 1959+650 
and Mrk421. On the other hand for non-orphan flaring 
this value is small i.e. Ay < 1 . 

In the flaring state, the proton luminosity L p for the 
highest observed proton energy E p = 168 TeV has to 
be smaller than the LEdd ~ 2.5 x 10 i6 erg s^ 1 and this 
gives T pl > 0.02 in the inner jet region. For our esti¬ 
mates we consider the hidden jet size R'f — 3 x 10 l0 
cm which is between R s and the blob radius R' b . This 
value of R'f corresponds to a day scale variability. We 
take day scale variability due to the fact that the flaring 
lasted for more than two weeks. If we consider small 
R'f then the seed photon density will increase as can 
be seen from ©. But this increase in density will not 
affect our results because we are using the scaling be¬ 
havior of Eq. ©. The constraint on the optical depth 
gives the lower limit on the seed photon density in the 
inner region n 1 f > 1.3 x 10 10 cm -3 . Again by assum¬ 
ing t' p y < t' d we obtain the upper limit on the optical 
depth T P y < 2 and this corresponds to the photon den¬ 
sity n' f < 1.3 x 10 12 c?n -3 . We have also estimated the 
seed photon density from Eq.®, for e 7 = 0.35 MeV, 
which gives n'y f < 8.9 x 10 lo cm -3 . The upshot of this 
analysis is that we get the constraint 1.3 x 10 10 cm ' 3 < 
n'y f < 8.9 x 10 lo cm~ 3 , which shows that the photon 
density in this region is high. This range of photon 
density corresponds to the optical depth in the range 
0.02 < T P y < 0.13 and the proton flux at E p = 168 TeV 
is below the FEdd ~ 1-24 x 10 -8 erg cm~ 2 s _1 . Due to 
the adabatic expansion of the inner blob, the photon den¬ 
sity will be reduced to n 7 and the energy will dissipate 
once these photons are in the bigger cone. So even if we 
have two-zones (the inner and the outer), only the outer 
zone will be responsible for the observed synchrotron and 
IC peaks. From the leptonic model fit to the SED, the 
magnetic field in the outer jet region is B' = 0.26 G and 
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higher magnetic field is expected in the inner jet region. 
The maximum proton energy in the hidden jet region 
will be E PjTnax ~ 10 18 (£>j/lG) eV and for larger mag¬ 
netic field the E PtTnax can even be higher. 

In the normal jet scanario we estimate the proton 
flux needed to explain the observed TeV 7 -rays. Cor¬ 
responding to the background energy e 7 = 0.35 the pho¬ 
ton density is n 7 ~ 4 x 10 3 cm ^ 3 and the the optical 
depth is r P7 ~ 1.4 x 10 -8 . As discussed above, for back¬ 
ground SSC photon energy 0.35 MeV, the observed TeV 
photon has the energy T 7 = 16.85 TeV and its flux is 
F 1 ~ 3.16 x 10~ n erg s~ cm~ 2 . The high energy proton 
responsible for the photohadronic process to produve this 
TeV 7 -ray has energy E p = 168 TeV and we can estimate 
its flux from Eq. m which gives F p ~ 10 6 x Fem ■ This 
shows that the normal jet model needs super Edding¬ 
ton power in protons to explain the high energy peaks, 
whereas, the inner jet scenario exterminates this extreme 
energy requirement. 

The high energy protons will be accompanied by high 
energy electrons in the same energy range. In the mag¬ 
netic field of the jet, these electrons will emit synchrotron 
photons in the energy range 4 x 10 19 Hz to 2 x 10 23 Hz , 
which is in the lower part of the SSC spectrum and will 
not be observable because of the lower flux. Also the SSC 
emission will take place from these high energy electrons 
and the energy of these IC photons will be Eic ~ lle syn . 
By considering the electron Lorentz factor in the range 
7 x 10 2 < 7 e < 4 x 1O 4 0 and the peak energy of the 
synchrotron photon e syn ~ 10 18 Hz, the SSC process can 
contribute in the energy range 2 GeV < Eic < 6.6 TeV. 
But the details of the SSC flux depends on the breaks 
in SSC spectrum and the spectral index. It is observed 
that during the flaring of Mrk 421, the X-ray emission 
reached the peak, days after the TeV emission, which 
poses a serious challange to the SSC model[llj]. The 
SSC contribution to the multi-TeV band will be very 
much supressed. Similarly the multi-TeV photons in the 
energy range 0.25 TeV < T 7 < 16.85 TeV can inter¬ 
act with the background photons to produce e + e _ pair 
and these individual electron or positron will have en¬ 
ergy E 1 / 2. To produce e + e _ pair the required thresh¬ 
old seed photon energy e 7 > 2 m 2 /E 1 is needed. Dur¬ 
ing the flaring, the multi-TeV 7 -rays in the energy range 
0.25 TeV < Ej < 16.85 TeV will interact with the soft 
seed photons in the energy range 0.05 eV < e 7 < 3.5 eV 
(in between the infrared and the visible range), where 
tr 77 ~ 1.7 x 10 ~ 25 cm 2 is the maximum cross section and 
for higher e 7 , the cr 77 will be smaller. The origin of these 
soft photons are from the synchrotron emission of 1 — 10 
GeV electrons in a magnetic held ~ 1 Gauss towards the 
base of the evolving jet as well as the ambient photons 
from the disk. On the other hand multi-TeV 7 -rays are 
produced beyond this region where the photons are in 
the low energy tail (0.35 MeV — 23.5 MeV range) of the 
IC photons. These two regions are distinct and there will 
not be enough low energy seed photons (0.05-3.5 eV) in 
the IC tail region. So, mostly the TeV photons will en¬ 


counter the IC photons and the pair production cross sec¬ 
tion for e 7 > 0.35 MeV is very small er 77 < 10 _3 ° cm 2 , 
which corresponds to amean free path A 77 > 10 19 cm. 
Hence, TeV photons will not be attenuated much due to 
e + e _ pair production. Also it has been observed that, 
during the flaring of Mrk 421, the variation in the light 
curves at optical and radio wavelengths are slight fllj. 
which shows that the low energy photon production was 
suppressed. The positron produced from the tt + decay 
will have energy E 1 / 2 and it will radiate synchrotron 
photons in the energy range 2 x 10 17 Hz to 9 x 10 2 ° Hz. 
The photon flux F e + syn from the synchrotron radiation 
of e + will be much smaller than F 1 (E 1 = 0.25TeV)/8, 
i.e., F e + syn <C 8 x 10 -11 erg cm~ 2 s _1 . This flux is well 
below the observed flux limit in the normal case as can be 
seen from Fig. [2] From the above analysis our conclusion 
is that, the photon fluxes from the synchrotron emission 
of the electrons and positrons are not observable during 
the flaring event of Mrk 421 in April 2004, which makes 
this flare orphan, like the one observed in 1ES 1959+650. 
In principle, the multi-TeV 7 -rays from the extragalactic 
sources can be reduced due to the absorption of TeV pho¬ 
tons by the diffuse extragalactic background light (EBL) 
through the process 'YtbV + 7 h —> e + e~ due to the energy 
dependent optical depth j35B37j . But for low redshifts 
and the energy rage of our interest, the optical depth 
does not vary much. Hence we can assume almost a con¬ 
stant optical depth[29j so that the spectral shape remains 
nearly unchanged. 

In conclusion, from the study of the flaring events in 
blazars, we deduce that the flaring phenomena can be ex¬ 
plained through the photohadronic interaction in a com¬ 
pact and confined region within the blazar jet where the 
photon density is high. From the study of 1ES 1959+650 
and Mrk 421, we deduce that orphan flaring can only be 
possible for those blazars which have a deep valley in be¬ 
tween the end of the synchrotron SED and the beginning 
of the SSC SED as shown in Fig. [2] We note that the 
HBLs Mrk 501 and PG 1553+113 are possible candidates 
for orphan flaring in future. For the orphan flaring events 
we find A 1 1 and for non-orphan flaring A 7 < 1 . 

If Mrk 501 were to produce an orphan flare then the 
flare energy will lie in the range ITeV < E 1 < 8.6 TeV 
and this corresponds to the background SSC photon en¬ 
ergy in the range 4.3 MeV > e 7 > 0.5 MeV. We esti¬ 
mate this by taking the parameters of Mrk 501 as fol¬ 
lows: V = 12, z = 0.034 from Ref.[3g]. In this case 
the maximum Fermi-accelerated proton energy will be 
E p < 10T 7 ~ 86 TeV. Also, another condition which 
must hold for the orphan flaring is 


T(8.6TeV) 

7 " <!>ssc(0.5 MeV) 



-a+3 

e 8.6/(E c /TeV) 


> 1 . 


(17) 

It has to noted that, if V changes, accordingly the values 
of Ej, e 7 and the E p will also change, but the condition 
of Ea. (fT71) is independent of these changes. 
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V. CONCLUSIONS 

The orphan flaring of Mrk 421 can be explained well 
by the hadronic model. We observe that, in this model, 
the multi-TeV photon flux is proportional to 'hssc(£ 7)1 
E~ a+3 and an exponential decay term as shown in Eq. 
(fT4l) . This implies that the Fermi accelerated protons in¬ 
teract with the background photons (in the low energy 
tail) of the SSC spectrum. During the April 2004 flar¬ 
ing of Mrk 421, we have shown that the flux from the 
synchrotron emission from the high energy e + and e~ is 


suppressed relative to the normal flux implying that the 
flaring was orphan in nature, like the one observed in 
1ES 1959+650. We have also shown what type of blazar 
spectrum will result in orphan flaring and predict that 
Mrk 501 and PG 1553+113 are possible candidates for 
orphan flaring in the future. 

Monitoring of these objects by the TeV gamma-ray 
telescopes will shed more light on the details on the or¬ 
phan flaring mechanism. 

The work of S.S. is partially supported by DGAPA- 
UNAM (Mexico) Project No. IN110815. 
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